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Elastic scattering of electron vortex beams on magnetic materials leads to a weak magnetic con¬ 
trast due to Zeeman interaction of orbital angular momentum of the beam with magnetic helds 
in the sample. The magnetic signal manifests itself as a redistribution of intensity in diffraction 
patterns due to a change of sign of the orbital angular moment. While in the atomic resolution 
regime the magnetic signal is most likely under the detection limits of present transmission electron 
microscopes, for electron probes with high orbital angular momenta, and correspondingly larger 
spatial extent, its detection is predicted to be feasible. 


Rapid developments in nanoengineering call for char¬ 
acterization methods capable to reach high spatial resolu¬ 
tion. In this domain, the scanning transmission electron 
microscope (STEM) provides a broad scale of measure¬ 
ment techniques ranging from Z-contrast [1] or electron 
energy-loss elemental mapping [ 2 ] , differential phase con¬ 
trast (DPC) [3, 4], via local electronic structure studies of 
single atoms [5] to counting individual atoms in nanopar¬ 
ticles [ 6 ]. As a specific case of high-spatial resolution 
electron energy loss spectroscopy, an electron magnetic 
circular dichroism (EMCD) method has been introduced 
[7] as an analogue to x-ray magnetic circular dichroism, 
which is a well established quantitative method of mea¬ 
suring spin and orbital magnetic moments in an element- 
selective manner [ 8 , 9]. 

Recenly, the introduction of electron vortex beams 
(EVB) [10-12], i.e., beams with nonzero orbital angular 
momentum, aimed at probing EMCD at atomic spatial 
resolution. It was shown theoretically that EVBs need 
to be of atomic size in order to be efficient for magnetic 
studies [13-15]. Several methods of generating atomic 
size electron vortex beams have been proposed [16-19], 
yet an experimental demonstration of atomic resolution 
EMCD has not been presented in the literature. 

An alternative route to utilizing EVBs for magnetic 
measurements is based on Zeeman interaction between 
their angular momentum and the magnetic field in the 
sample. The Pauli equation for an electron with energy 
E in an electrostatic potential V (r) and a constant mag¬ 
netic field Be reads 


tween the magnetic field and the orbital and spin angu¬ 
lar momenta of the electron beam. A previous study has 
indicated that the effect of spin on elastic scattering is 
very weak [20]. Moreover, generating intense spin polar¬ 
ized electron beams remains a technological challenge [ 21 ] 
and so far magnetic field mapping with spin-polarized 
electrons in the TEM could not be demonstrated. While 
the spin angular momentum of electrons in the propaga¬ 
tion direction is at most |, EVBs can be generated with 
very high orbital angular momenta (0AM) [ 12 , 22 , 23], 
which permits an increase of the Zeeman interaction by 
more than two orders of magnitude. 

In this Letter, we show that there is a magnetic con¬ 
trast in elastic scattering of EVBs originating from the 
enhanced Zeeman interaction of the beam 0AM with 
magnetic fields in the sample. The described effect is 
sensitive to magnetic fields parallel to beam-direction, 
which would complement holographic or DPC methods 
measuring the in-plane components of the magnetic field. 

Eor a realistic description of magnetism in a solid, tak¬ 
ing into account merely a constant magnetic field is insuf¬ 
ficient. Hence, we consider a stationary Pauli equation 
with a non-uniform magnetic field [24] B(r) = V x A(r) 
and corresponding vector potential A(r) in Coulomb 
gauge, V • A(r) = 0. Due to large acceleration voltages 
commonly applied in TEM, a relativistically corrected 
electron mass m = qmo is used. Subsequently, as in the 
derivation of the conventional multislice method [25] , we 
introduce a paraxial approximation [26] via the substitu¬ 
tion 
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where —e is the electron charge, m is the electron mass, 
p = —ihV is the momentum operator, L and S are the or¬ 
bital and spin angular momentum operators, and 4^(r) is 
a two-component spinor wavefunction. The second term 
on the left hand side of Eq. 1 manifests a coupling be¬ 


and neglect the second derivatives of the envelope func¬ 
tions - 0^4 (r) with respect to the beam propagation direc¬ 
tion k = (0,0, k). The resulting two-component paraxial 
Pauli equation reads [27] 
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which upon setting A = B = 0 reduces to the parax¬ 
ial Schrodinger equation [26] for each of the spin com¬ 
ponents separately. Eq. 3, however, represents a 
system of two differential equations coupled via an inter¬ 
action of the spin of the probe with the magnetic field in 
the sample. It can be integrated slice-by-slice according 
to [25, 28] 



'0(x, y,z-\- Az) = Z{eJ^z ~ 

^ (4) 
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where Z is Dyson’s 2 ;-ordering operator. Similar com¬ 
putational methods were recently discussed for the fully 
relativistic case of the Dirac equation [20] and in the con¬ 
text of spin polarization devices [29]. 

In the following numerical simulations we constructed 
the electrostatic potential from tabulated values of in¬ 
dependent atoms [26], whereas the magnetic vector po¬ 
tential and the corresponding magnetic field are obtained 
from density function theory (DFT) in the following way. 

In a crystal, B consists of a constant part due to the 
saturation magnetization Be = /ioMg and a periodic part 
Bp that averages to zero. The constant part originates 
from a non-periodic component of vector potential Anp = 
X r, while the periodic part of the magnetic field 
originates from Ap computed as a periodic solution of 
AAp(r) = —where j(r) is the spin current density. 
By using a Gordon decomposition and neglecting orbital 
currents, the spin current density is [20] j(r) = V x m(r), 
where the spin magnetization density m(r) is computed 
from electronic structure spin DFT calculations. 

We expect that the microscopic variations of the mag¬ 
netic field. Bp, will only play a role in atomic resolution 
regime. For larger probes, such as EVBs with high 0AM, 
effects of these variations average to zero and only the 
constant part of the magnetic field Be will influence the 
scattering on top of the Coulomb potential. The situa¬ 
tion can then essentially be understood in terms of Eq. 1. 

Results of the procedure described above, applied to 
bcc Fe, are illustrated in Fig. 1. In the case of collinear 
magnetism, m(r) is parallel to the 2 ;-direction, whereby 
j(r), and in the gauge chosen here also A(r), have non¬ 
zero X- and ^-components only. The spin magnetiza¬ 
tion density as obtained via a collinearly spin-polarized 
full-potential linearized augmented plane wave [30] cal¬ 
culation in the generalized gradient approximation [31] 
is shown on an x^-cross section containing the central Fe 
atom of the bcc unit cell. Fig. la. The ^-component of 
the B-field, Fig. lb, and the x-component of the Ap-field, 
Fig. Ic, are plotted within the same plane. Note that the 



Figure 1. a) Spin magnetization density, b) ^-component of 
the B-field, c) x-component of the Ap-field and d) the B-field 
in a unit cell of bcc Fe as obtained from DFT calculations. 

^-component of Bp reaches values in the order of 60 T— 
significantly larger than jHoMs = 2.2 T in bcc Fe. Ay (not 
shown) is identical to A^ rotated by 90° about the z-axis 
and Az = 0 everywhere. Finally, the microscopic Bp- 
field (to which a 5c = 2.2 T in the ^-direction should be 
added) is plotted as a vector field in one unit cell. Fig. Id. 
Although the shape of the spin density is very similar to 
that of Bz, they are not identical, and, even though only 
collinear spin density along the ^-direction is considered, 
the Bp-field has non-zero x- and ^-components. 

Due to persisting limitations in the creation of EVBs, 
electron beams with large OAMs (;» Ih) cannot be fo¬ 
cussed on an atomic scale. We thus concentrate first on 
a situation, where we do not aim for atomic resolution, 
but rather enhance the Zeeman interaction by a large ini¬ 
tial 0AM (= Ih) of the beam. Using Eq. 4 we propagate 
electron beams with an initial 0AM of 20fi, 30h and 40fi, 
through a bcc Fe crystal of thickness up to 400 unit cells 
(115 nm). The radial shape of the beams is described by 

i’i{k±,(p) e'’-'^Q{ak - k±), (5) 

where k± and 0 are cylindrical coordinates in k-space and 
a is the convergence semi-angle. The lateral supercell 
dimension was 48 x 48 unit cells and each unit cell was 
discretized on a 64 x 64 x 64 pixel grid. The acceleration 
voltage was 200 kV and a = 10 mrad, corresponding to 
outer full-widths at half-maximum of 2.2, 3.1, and 3.9 nm 
for I = 20, 30, and 40, respectively. Beams were centered 
on an atomic column, but we have verified that the results 
do not depend on the exact beam position, as is expected 
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for beams with spatial extent significantly larger than the 
crystal unit cell. 



Figure 2. a) Magnetic signal (difference of intensities) as a 
function of collection angle. Differences are shown for unpo¬ 
larized beams with I = ±20, ±30 and ±40, as well as beams 
with I fixed to —20 or 30, but with difference taken over oppo¬ 
site spin channels. For I = ±20 the total signal (sum of inten¬ 
sities for opposite 0AM) is also shown. Sample thickness was 
set to t = 42 unit cells (12 nm). b) Thickness dependence of 
the relative magnetic signal for a circular collection aperture 
with collection angle indicated (/ = ±30). 

A non-spin-polarized beam is in a mixed state, with 
50% of electrons with spin-np and 50% spin-down, re¬ 
spectively. Therefore each simulation consists of two 
runs, one for each spin orientation, and the resulting 
diffraction patterns were averaged over the two spin ori¬ 
entations. It is worth mentioning that the proportion of 
spin-up electrons scattering into spin-down states or vice 
versa is negligible (of the order 10“^^), although it has 
been suggested that for magnetization in the x^-plane 
the spin-flip scattering can be more significant [32]. 

The Zeeman interaction leads to a redistribution of 
intensity in the diffraction pattern. Total intensity of 
scattered electrons is of course the same for positive or 
negative 0AM, but the intensity of electrons scattered 
to smaller or larger angles varies depending on the 0AM 
(see Fig. 2a). For large collection angles the intensity of 
scattered electrons saturates at a value of one, as dic¬ 
tated by normalization of the initial probe wavefunction. 
The intensity difference for, e.g., I = ±20 shows a peak 
around a collection angle of 10 mrad, after which its am¬ 


plitude decreases eventually reaching zero. Computing 
such differences in a simulation with zero magnetic fields 
merely yields a numerical noise around ten orders of mag¬ 
nitude smaller. The kink observed close to ^ = 130 mrad 
appears near a higher order Lane zone [33]. 

Notice, how the magnetic signal is approximately pro¬ 
portional to the initial 0AM of the EVB. This suggests 
that the strength of this signal can be further scaled up 
for beams with larger 0AM. Consistently with Eq. 1, the 
magnetic signal obtained as an intensity difference from 
opposite spin channels, but for a fixed initial 0AM, is 1) 
independent of 0AM, and 2) of the same order of magni¬ 
tude as the intensity difference due to changing the sign 
of 0AM, when normalized per unit of 0AM. 

In Eig. 2b) the magnetic signals, for a disc shaped re¬ 
gion with collection angles indicated in the legend, are 
shown for the I = ±30 case as a function of sample thick¬ 
ness. Accordingly, the magnetic signal can become sig¬ 
nificantly stronger for thicker samples. After 100 u.c. 
significant values of the relative magnetic signal in the 
range of 10“^ are observed. Considering that the sig¬ 
nal is proportional to 0AM and OAMs of size several 
hundreds have been reported [23] , signal strengths of few 
percent can be reached. While the intensity differences 
sensitively depend on the shape of the probe, magnetism 
should be measurable in an experimental setup with sym¬ 
metrical ±/ and —I 0AM probes — such as those gen¬ 
erated by holographic zone plates [11, 12]. Relative sig¬ 
nals of several percent are well within the the detection 
limits of current bright field EVB STEM experiments, 
which easily integrate order of 10^ electrons (for typi¬ 
cal probe currents, dwell times, collection angles) with a 
corresponding shot noise order of 10^, i.e., l%o. 

Now we turn our attention to the atomic resolution 
regime. As was mentioned above, to focus EVB with 
large 0AM onto atomic scales requires very large con¬ 
vergence angles in the range of hundreds of milliradi- 
ans, which is outside of present instrumental possibili¬ 
ties. Therefore, we restrict ourselves to EVBs of small 
initial 0AM, thereby reducing the magnitude of the Zee- 
man term. However, we remind that the strength of local 
microscopic magnetic fields can reach substantial values 
of several tens of Teslas (see Eig. lb), which could po¬ 
tentially lead to significant magnetic signals even in the 
atomic resolution regime. To assess the effect, we have 
performed simulations using beams with / = ±1 and a 
rather large convergence angle of 40 mrad at an accel¬ 
eration voltage of 300 kV. The supercell dimension was 
24 X 24 X 100 unit cells of bee Ee (28.7 nm thick), each 
cell discretized on a 112 x 112 x 112 grid. The collection 
angle was set to 5 mrad. 

Eig. 3 summarizes atomic resolution simulations. The 
STEM image (Eig. 3a) should be compared to the 
magnetic signal computed from the / = ±1 difference 
(Eig. 3b) and spin-difference at / = 1 (Eig. 3c). Note 
that a mirror image of vortex beam with 0AM equal to 
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Figure 3. Simulated STEM image for a collection angle 
5 mrad, as obtained from a) the spin averaged beam for 
I = +1, b) magnetic signal obtained as difference signal due 
to 0AM I = ±1, c) difference between spin-up and spin-down 
beams for / = +1. 


Ihis a vortex beam with 0AM equal to —Ih [34]. For this 
reason the magnetic signal is obtained as a difference of 
intensity of / = +1 beam at position {x^y) and / = — 1 
beam at (^, x). 

Both spin and 0AM differences are of the same order 
of magnitude, which is about 10“^ of the total signal. In 
the following we explore some routes to enhance the mag¬ 
netic signal. An increase in its relative strength can be 
achieved by a further optimization of sample thickness, 
convergence angle, acceleration voltage, initial 0AM and 
collection angle. To find the global maximum of this 
multidimensional optimization problem is a formidable 
task, which will furthermore depend on the material to 
be investigated. Therefore, we focus on a few parame¬ 
ters only and assess the increase in magnetic signal in an 
approximative manner. 

To investigate the effect of larger 0AM on the mag¬ 
netic signal strength, a few beam positions were recalcu¬ 
lated with with I = ±2 as well as I = ±4, see Fig. 4a)-f). 
Accordingly, the proportionality between magnetic signal 
and I is lost in the atomic resolution regime. This is most 
likely due to a different form of the magnetic interaction 
oc A • p, which is not anymore directly proportional to 
the 0AM and depends on details of spatial distribution 
of the magnetic field and the probe wavefunction. Note 
for example the radial intensity profiles for beams with 
I = 1, 2,4. The differences are mostly due to strong pin¬ 
ning of beams with low 0AM to atomic columns [35], 
less pronounced for 0AM=2/i and 4/i, respectively. Nev¬ 
ertheless, magnetic signals are still somewhat stronger 
for beams with larger 0AM. 

An alternative route to enhancing the magnetic signal 
consists of reducing the acceleration voltage. Upon in¬ 
spection of Eq. 3 we note an additional prefactor 7 “^ in 
front of the magnetic coupling compared to the electric 
one resulting in a relative increase of the magnetic signal 
at lower acceleration voltages. Fig. 4g)-i) compares re¬ 
sults obtained with voltages 100 kV and 300 kV for beams 
with initial 0AM of and other parameters kept fixed. 
An increase of magnetic signal by a factor of ^ 3 or 4 
can be observed for the lower acceleration voltage. 



0(mrad) 


Figure 4. Sum in a)-c) and difference in d)-f) of signal inte¬ 
grated over collection angles from 0 to ^ for +/ and —I beams 
with / = 1, 2 and 4 at beam positions (0,0) in a), d), g), 
(1,1)^ in b), e), h) and (1,0)f in c), f), i). g)-h) show 

the difference between acceleration voltage of 100 keV or 300 
keV with other parameters kept same for the three respective 
beam positions. 


Combining all the effects a further optimization of all 
parameters is suited to increase the relative magnetic sig¬ 
nal strength by one order of magnitude to 10“^. Yet, 
the relative magnetic signal strength of 10 “^ means that 
it will be extremely sensitive to scan noise, drifts and 
changes of sample orientation during data acquisition, 
which renders atomic resolution measurements of a mag¬ 
netic signal based on the Zeeman interaction of 0AM 
with magnetic fields in the sample extremely challeng¬ 
ing, most likely beyond the possibilities of present in¬ 
struments. 

In conclusion, we have demonstrated computationally 
that the elastic scattering of electron vortex beams on 
magnetic samples in the TEM depends on the relative ori¬ 
entation of the initial 0AM and the magnetization in the 
sample. In principle, this effect opens a new way for mea¬ 
surement of magnetic properties. For beams with 0AM 
of few hundreds of the predicted relative strength of 
magnetic signal should reach up to a few percent, call¬ 
ing for an experimental verification. If successful, this 
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permits a new way of characterization of magnetic prop¬ 
erties at about 10 nm spatial resolution. In the atomic 
resolution regime, the calculated relative magnetic signal 
strength reaches only up to 10“^, making it unlikely to 
be detected with present-date instruments, particularly 
due to scan noise and unavoidable sample drifts. 
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